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Abstract. We tested the ejector-propeller model of the Be/X-ray binary LSI-l-6l''303 by using the parameters predicted by 
the model in the calculations of the X-ray and radio variability. The results are: (1) in terms of the Ejector-Propeller model, the 
X-ray maximum is due to the periastron passage; (2) the radio outburst can be really a result of the transition from the propeller 
to the ejector regimes; (3) the radio outburst will delay with respect to the X-ray maximum every orbital period. The proposed 
scenario seems to be in good agreement with the observations. 
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1. Flaring beliaviourof LS 1+61 °303 

' LS lH-61°303 (V615 Cas, GT0236H-620) is a well known Be/X- 
ray binary with periodic radio outbursts every 26.5 d, assumed 

' to be the orbital period. After the discovery of such periodic 
events (Gregory & Taylor 1978; Taylor & Gregory 1982) and 
the first model (Maraschi & Treves, 1981), this massive sys- 
tem has been studied extensively since nearly two decades 
ago. Hereafter, we will use the late radio phase ephemeris 
P = 26.4917 d and phase zero at JD2443366.775 (Gregory, 
Peracaula & Taylor, 1999). 

The flaring radio emission of LS 1h-61°303 has been mod- 
eled by Paredes et al. (1991) as synchrotron radiation from an 
expanding plasmon containing relativistic particles and mag- 
netic fields. The dependence of the radio outburst flux density 
on the frequency, the time delay in the outburst peak, and the 
general shape of the radio light curves are in general well ac- 
counted for by continuous injection of synchrotron emitting 
electrons into the plasmon volume. However, the formation de- 
tails of such plasmon were beyond the scope of the Paredes et 
al. (1991) model. In fact, its initial physical properties were 
considered as the free parameters to be fitted. According to 
Zamanov (1995), the genesis of such plasmon can be inter- 
preted as a result of a transition in the accretion regime of the 
neutron star as it probes different parts of the Be envelope in 
an eccentric orbit. The proposed transition would be from so 
called propeller (P) to the Ejector (E) regime. In this context, 
propeller means accretion onto the magnetosphere and ejector 
is often named "young radio pulsar". 



Our purpose here is to try to check the agreement between 
the plasmon parameters predicted by the theory of a PE tran- 
sition and those required to reproduce the observed radio light 
curves using the Paredes et al. (1991) model. We advance here 
that the agreement seems to be good at least qualitatively, in- 
cluding the possibility to explain the X-ray luminosity and vari- 
ability. This agreement provides a consistent physical interpre- 
tation for the outburst mechanism in the LS 1h-61°303 system. 

2. Simultaneous Radio and X-ray observations 

The variability of the source in the X-ray domain was detected 
in ROSAT observations (Goldoni & Mereghetti, 1995). On two 
occasions, the LS 1h-61°303 flaring events have been moni- 
tored by ground based radio facilities and space X-ray satellites 
in a coordinated way. The results were presented in Taylor et 
al.(1996) and Harrison et al.(2000). The behavior observed is 
reproduced in Fig. |l], using the latest ephemeris determination. 
The observations during 1992 are obtained with ROSAT (0.5-2 
keV) and die VLA (1.5 GHz and 4.9 GHz), and those during 
1996 - with RXTE (2-10 keV) and GBI [] (2.25 GHz and 8.3 
GHz). From this figure, it is clear that in both cases the radio 
peak is delayed with respect to X-ray maximum. 

3. Mass accretion rate 



' The Green Bank Interferometer is a facility of the USA National 
Science Foundation operated by NRAO in support of the NASA High 
Energy Astrophysics programs. 
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Fig. 1. The simultaneous X-ray and radio observations: 

(a) and (b) - August 1992, Taylor et al. (1996) 

(c) and (d) - March 1996, Harrison et al. (2000) 

A clear phase shift between X-ray and radio is visible in both 

cases. 

In (b) the triangles refer to 1.5 GHz and the crosses - to 4.9 
GHz flux density. The X-ray flux is everywhere in units 10^^^ 
erg cm~^ s^^. 

3. 1. Bondi-Hoyle accretion 

Usually, to estimate the mass accretion rate onto a compact ob- 
ject in wind fed binaries the Bondi-Hoyle-Littleton formula is 
used in the form: 

Ma = TT —3 p, (i) 

^rel 

where M^s is the mass of the neutron star, Vrei is the NS ve- 
locity relative to the surrounding matter, p is the density in the 
surrounding matter 

The structure of the Be star disks are not yet well under- 
stood, even for the most important parameters, such as the ve- 



locity and density laws, different assumptions have been pro- 
posed. Following Waters et al. (1988), we will accept a density 
structure in the form 

p = Poir/R,)-'' and Vr ^ Voir/ R,)"-^ , (2) 

where R^ is the radius of the Be star and Vr is the outward 
velocity. The density parameter n can take values n = 2 — 4 
for different objects. The initial velocity, Vq (the velocity at 
r = R^,), most probably lies between 2 and 20 km s^^. Finafly, 
a Keplerian law is adopted concerning the disk rotation: 

V^r) = ^GMi/r (3) 

where Mi is the mass of the primary component in the system, 
i.e. the Be star 

As a result of the motion of the neutron star throw the 
Be circumstellar disk the mass accretion rate will vary along 
the orbital period. A few examples, calculated using different 
parameters for the Be circumstellar disk, are shown in Fig.|[ 
During the calculations we assumed always an orbital eccen- 
tricity e — 0.6, a neutron star mass of M^s = 1.4 M©, 
Ml — 10 Mq and an orbital period Porb=^6-^9n days. 

As it is visible in Fig.|[ we can obtain very strange mass ac- 
cretion rates in some cases. In Fig.^(b,c), the average mass ac- 
cretion rate is considerably higher than the mass loss rate of the 
Be star. This is a false result of the used formula. The Bondi- 
Hoyle description is valid for accretion from infinite medium. 
It assumes that all the matter entering into a cylinder with ra- 
dius equal to the accretion radius (Ra — \/'2GMns /Vrei) ^^^^ 
be accreted. In the context of a binary star, this is only valid 
provided that 

- the motion is high supersonic V^re; >> Cg, 

- Ra « r, 

- the whole accretion cylinder is filled with the outflowing 
wind. 

Unfortunately these conditions are not always fulfilled in 
Be/X-ray binary systems. 

3.2. Mass accretion rate in Be/X-ray binaries 

As it was shown in the previous section the Bondi-Hoyle for- 
mulae are not giving good results in case of Be/X-ray binaries. 
The exact mass accretion rate could be found only if we con- 
sider simultaneously the formation of the wind and the influ- 
ence of the compact object onto the wind. Here we propose a 
modification of the mass accretion rate formula that is appro- 
priate for Be/X-ray binaries. In order to do so, the Bondi-Hoyle 
formula has to be rewritten in a way such that the accretion ra- 
dius explicitly appears: 

Ma = 7rRlpVrel, (4) 

where the density in an outflowing wind, with disk opening 
angle 9, is given from the continuity equation: 

Mioss = 47rr2 sin 9 p, (5) 

The wind in the disk is most probably subsonic (i.e. 
Okazaki, 2001). A keplerian rotation in the disk will give in 
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Fig. 2. Mass accretion rate using the Bondi-Hoyle formulae. Everywhere the orbital parameters are Mi — 10 Mq, Mjvs = 
1.4 Mq, Porb=^^-^9ll d, e — 0.6. The wind parameters used and the calculated average mass accretion rate are : 

(Al) Vo=5 km s"\ n = 3.25, non-rotating disk, Ma = O.ll Mioss 

(Bl) Vo=5 km s~\ n = 3.25, keplerian disk, M„ = 2.9 Mioss 

(CI) Vo=10km s-\ n = 2.1, keplerian disk. Ma 9.7 Mioss 

The use of the Bondi-Hoyle formula gives obviously false values( at least in Al and Bl ) values of the average mass accretion 
rate: 

(A2), (B2), (C2) - the same parameters, but the factor of the filling of the accretion cylinder included. 

(A3), (B3), (C3) - the same parameters, but the upper limit set Ra < r/A, and the filling factor included. The corresponding 
average mass accretion rates are Ma/ Mioss — 0.09, 0.08, 0.36, for A3, B3, and C3 respectively. 



many cases a low relative velocity between the NS and the 
wind. This means that the accretion radius depends also on the 
speed of sound, Cs : 



Ra — 



2GMns 
VI, + c2 



(6) 



Using this formula (with or without term), the accre- 
tion radius of a neutron star can achieve values bigger than the 
distance between components. However the disks in the Be/X- 
ray binaries are more or less similar to the disks in the "normal" 
Be stars (Zamanov et al. 2001) although some differences exist. 
This means that the accretion radius can not be very big, oth- 
erwise the NS would control the motion of the wind and this is 
impossible because the mass ratio is Ml /Mats > 4 (typically 



Ml > 8 Mq and M, 



1.4Mo) We will thus adopt here 



an upper limit Ra 
Be star 



< 0.25 r, where r is the distance from the 



Because the outflowing wind is confined in a disk, it is pos- 
sible that the vertical size of the disk is less than the accre- 
tion radius. Consequently, the accretion cylinder would not en- 
tirely filled by the wind material. In case that Ra > y where 
y is the vertical size of the outflowing disk y — rsinO, the 
accretion rate must be reduced with pure geometrical factor 
2^R-^[y./W^ + Rlsin-\y/Ra)]. 

We believe that the proposed modifications give consider- 
ably better results concerning the mass accretion rate, i.e. the 
NS will accrete a fraction of the wind. The results of the calcu- 
lations are demonstrated in Fig.|| A3,B3,C3. 
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Fig. 3. A sketch of the Ejector-propeller model. The dashed lines indicates the orbit of the neutron star. Upper - accretion onto 
the magnetosphere at periastron, down - ejector and the cavern around the neutron star 



4. X-ray variability in thie iighit of 
Ejector-Propeiier modei 

A sketch of the ejector-propeller model of LS Ih-61°303 is il- 
lustrated in Fig.^. Following the suppositions the neutron star 
will act as propeller (accretion onto the magnetosphere) at peri- 
astron and as ejector (young radio pulsar) at the apastron. The 
transitions between the regimes will depend how material is 
captured in the the neutron star accretion cylinder, possible for- 
mation of accretion disk during the accretion onto the magne- 
tosphere, and even the possible flip-flop instability of accretion 
(Livio etal 1991). 

Following the assumptions of the ejector-propeller model, 
a lower mass capture rate is expected during the ejector regime. 
In contrast, higher mass capture rates are expected during the 
the propeller regime close to periastron passage. During the 
propeller stage the X-ray luminosity can be estimated as 



LxiP) = 22/7(GA/)8/^M-4/7^y9/7^ 



(7) 



where G is the constant of gravitation, M is the mass of the 
neutron star, /i is its magnetic dipole moment (which is adopted 
jj, — 2 X 10'^°G cm^^), and Ma is the mass accretion rate. 

During the ejector stage efficiency of conversion of the 
spin-down energy into high energy emission can be expected 
to be about 1-10% (Tavani et al. 1996). The estimated spin- 
down luminosity of the ejector is expected about 1 — 5 x 10'^^ 
ergs^^, so as reasonable value can be adopted Lx{E) — 
2.1033-2.1034 ergs-i. 

The other possibility to estimate the X-ray luminosity of the 
expanding plasmon is to consider Bremsstrahlung emission of 
the shocked gas. In this case Taylor et al. (1996) also obtained 
2 X 10^* ergs^^. The both estimates are in good agreement 



and are similar to the expected values of the magnetospheric 
accretion. 

It deserves to be noted that in case of transition between 
accretor and propeller a considerable luminosity gap will exist, 
about 100 times difference in the X-ray luminosity, is expected 
and observed (Corbet 1996; Raguzova & Lipunov, 1998). But 
this is not the case at the transition between ejector and pro- 
peller, where the X-ray luminosities are of the same order So, 
we will assume that Lx (ejector) = Lx (propeller) not to 
complicate our qualitative results. 

The X-ray luminosity calculated in this way is plotted in 
Fig.^ The orbital phase zero corresponds to the periastron, the 
adopted mass loss rate of the Be star is Mioss = 10^^ Mq 
yr^^. The calculated X-ray variability is very similar to the ob- 
served (Fig. [l|). The X-ray maximum is at periastron and this 
means that the periastron passage corresponds to radio phase 
^ 0.5 using the modern radio ephemeris. 



5. Plasmon parameters and radio outbursts 

The switch on of the ejector will create a cavern around the neu- 
tron star. This cavern will start to expand under the pressure of 
the relativistic wind of the ejector and can be identified with the 
expanding plasmon, which radio emission successfully fit the 
observed radio outbursts of LS Ih-61°303 modeled by Paredes 
et al. (1991). The question here is: are the parameters expected 
at the propeller-ejector transition appropriate for the plasmon? 

The most important parameters which appeared in the ra- 
dio model of Paredes et al. (1991) are initial radius (Rq), ex- 
pansion velocity (Vc^p), initial magnetic field, injection time 
interval, power law index of electrons, and the injection rate. In 
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Fig. 4. The calculated X-ray luminosity of the the neutron star in LS 1+61 °303 using different wind parameters. The Lx is 
calculated as magnetospheric accretion luminosity and supposing that Lx (ejector) = Lx {propeller). The radio outbursts are 
expected to peak after the X-ray maximum, when the mass accretion rate is lower, and the NS can change its regime from 
propeller to ejector. 

(a) Vo=5 km s~ \ n = 3.25, non-rotating disk, Ma = 0.09 Mjoss 

(b) Vo=5 km s"\ n = 3.25, keplerian disk. Ma = 0.08 Mioss 

(c) Vo=10 km s"\ n = 2.1, keplerian. Ma = 0.36 Moss 
id) Vo=10 km s"\ n = 3.25, keplerian. Ma = 0.03 Moss 

(e) Vo=2 km s"\ n = 3.25, keplerian, Mg = 0.20 Moss 

(f) Vo=5 km s"\ n = 3.75, keplerian, Ma_= 0.02 Mioss 

(g) Vo=5 km n = 2.25, keplerian, Ma = 0.53 Moss 

(h) Vo=2 km s"\ n = 3.25, keplerian. Ma = 0.20 Mioss 

Everywhere we adopted mass loss rate in the outflowing disk Mioss = Mq yr~^ . 



the terms of Ejector-propeller model we can put constrains on 
the allowed ranges for this parameters: 

- 1. The time of plasmon appearance - the plasmon will 
appear after the periastron passage at low mass accretion rate 

(for more details see Zamanov, 1995 and references therein). 
So it means that the radio outbursts will peak always with some 



phase shift after the change of the regime from propeller to 
ejector. 

- 2. The initial radius expected at the propeller-ejector tran- 
sition will be of the order of light cyUnder radius, and in any 
case less than the accretion radius, which means that acceptable 
values are of the order of 0.01 — 3.0 Rq 
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- 3. The expansion velocity of the cavern in terms of 
the propeller-ejector transition can be estimated as (Zamanov, 
1995): 

V^xp = + \/L,nV^/{cMloss) (8) 

where c is the speed of light, Mioss and Vw are the mass loss 
rate and the wind velocity of the Be star respectively, Lm is 
magneto-dipole luminosity of the young radio pulsar Adopting 
14, = lOOkms^i for the disk wind and = 1000 km s^^ for 
the polar wind and L„j = 1 — 5 x lO'^^erg s^^ the cavern expan- 
sion will be with velocity 200 — 1300 km s^^. In the radio light 
curves we will adopt values V^oxp =200, 400, 1200 kms^^. 

- 4. Injection time interval - it is the time neutron star acts as 
ejector The switch on of the ejector will be when the neutron 
star goes away from the B star and the switch off will be when 
it enters in the dense circumstellar disk approaching the peri- 
astron. From the mass accretion rate behavior (Fig.^ right and 
Fig.Q) we can expect injection time interval 0.4 — 0.7 x Port or 
in other words 10 — 20 days. 

- 5. Injection rate - the total injection rate of the Crab into 
the nebula at a spin period of 0.033 s^^ is estimated to be 
about 10"'°-10'*^ particles on base of the high-energy spec- 
trum (De Jager & Harding, 1992). From one side the spin pe- 
riod of the neutron star in LS I+61°303 required from the E-P 
model is longer than Crab. We expect Pspin = 0.15 — 0.20 s 
(see Zamanov 1995 for more details). It means ~ 100 — 1000 
times lower rate of injected relativistic particles (L„j ~ P^j^n)- 
From other side as a result of the magnetosphere accretion 
there are probably more particles for injection. We will adopt 
5xl0-38e-s-i 

- 6. Magnetic field - the origin of magnetic field inside the 
plasmon will be something like to the magnetic field of the 
Crab nebula but some additional contribution from the B star 
is possible. The magnetic field referred in the "radio" figures 
(Fig.^ ) are very high because they refer to the magnetic field at 
the initial plasmon radius. When the plasmon expands to sizes 
of 1 AU, as measured with VLBI, the magnetic field is then 
close to 1 Gauss in agreement with VLBI estimates based on 
equipartition arguments. (Note: Changes in the radius and the 
magnetic field do no alter the light curves provided that B 

is constant - magnetic flux conservation). 

- 7. Electron power law index is adopted j3 ~ 1.6, in orderto 
be consistent with observed non-thermal (i.e. negative) spectral 
indices when the plasmon becomes optically thin. According 
to synchrotron theory, the optically thin spectral index is a = 
(1 — p)/2. In our case, we have then a — —0.3 in agreement 
with typical observations. 

- 8. Energy range of relativistic electrons: 
TTieC^ < E < 10~^ erg 

- 9. Distance from the Earth: 2 kpc 

Using the above parameters appropriate for the ejector- 
propeller model and the plasmon prescription we generated ra- 
dio light curves. They are shown in Figjs} The radio behavior 
(amplitude and shape) is very similar to the observed in the 
radio monitoring of LS I+61°303. 



6. Discussion 

Using the parameters appropriate to the ejector-propeller model 
we calculated the X-ray and radio variability of the Be/X-ray 
binary LS 1+61 °303. The general shape of the radio behavior 
(amplitude and shape) is very similar to the observed in the 
radio monitoring of LS I+61°303, as it was already shown in 
the "plasmon" calculations by Paredes et al.(1991). The radio 
outbursts peak can be achieved 2-8 days after the appearance of 
the plasmon (see Fig.|| ). The appearance of the plasmon can be 
expected when the neutron star emerges from the denser part 
of the disk, i.e. 2—5 days after the periastron (orbital period 
is 26.5 days). In the calculated X-ray curves the X-ray maxi- 
mum corresponds to the periastron. The switch on of the ejec- 
tor will be 0.1 — OAPorb later It means that the radio peak 
will delay with 3-13 days after the X-ray maximum. This is 
the behavior observed in both cases of the simultaneous radio 
and X-ray observations (Fig.|l]). As can be seen in Fig.^ the 
total X-ray luminosity expected in Ejector-Propeller model is 
of the order 1 — 200 x 10"^'^ erg s^^ This is in good agreement 
with the observed by ASCA (Paredes et al. 1997) variability 
1 - 6 X 10^4 erg s"! in the 2-10 keV band. 

6. 1. High resolution radio maps 

Recent high resolution radio maps have evidenced a one sided 
radio jet at milliarcsecond scales (Massi et al. 2001). These 
authors interpret it as a microquasar bipolar ejection with sig- 
nificant Doppler boosting effect. The ejector-propeller model, 
explored in this paper, provides however an alternative inter- 
pretation. Considering that the plasmon will be formed in one 
side of the Be star (the apastron vicinity), a one sided radio jet 
is naturally expected. 

6.2. Asymmetry oftlie cavern 

To model the radio light curves we assumed that the plasmon 
expansion is spherically symmetric. But in the real case this 
region will be not symmetric. Different forms of the cavern 
are possible - closed or open, for more details see Lipunov 
& Prokhorov (1984). In addition, the fact that the neutron star 
ejects relativistic particles at the apastron of an eccentric orbit 
will lead to a radio source which is elongated in the direction of 
apastron. In any case, the radio source formed around the eject- 
ing neutron star will be elongated in the direction controversial 
to the direction to the B star. 

6.3. Strong and weak radio outbursts 

It was discovered that radio outburst peak flux density varies 
over a time scale of 1600 days. On the same time scale the 
Ha emission of the outflowing disk (Zamanov & Marti, 2000) 
varies as well, and may be even the X-ray maximum (Apparao 
2001). Our preliminary tests evidenced that a slower expansion 
velocity and stronger magnetic field give higher radio peak flux 
densities and the outburst peaks later Also a faster expansion 
velocity and lower magnetic field will result into lower peak 
flux densities and the earlier outburst peaks. These effects and 
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their connection with the conditions in the Be disk need of care- 
ful modeUng. 

6.4. Multiple Outbursts 

A mechanism of multiple formation of caverns around an eject- 
ing neutron star is proposed (Lipunov & Prokhorov, 1983). We 
can speculate that, in case of multiple peaks in the LS 1+6 1 °303 
radio outbursts, we may be observing such a multiple cavern 
formation. Another possibility is a density structure in the wind 
of the Be star, i.e. rings with higher density that can change the 
regime more than once during an orbital period. 

6.5. Gamma rays 

During the Ejector stage the gamma ray emission is thought to 
originate in the shock front at the boundary of the pulsar and 
stellar winds (Kniffen et al. 1997). Gamma ray production is 
possible during the propeller action too. Possible mechanisms 
are double layer formation in the interblob plasma (Wynn, King 
& Home, 1997) or high voltage at the boundaries of vortex 
structures (Wang and Robertson, 1985). 

6.6. Other possibilities 

The fact that the radio outbursts are not exactly periodical point 
to they are a result of a process which has appropriate condi- 
tions to begin somewhere at radio phases 0.5-1.0, but when the 
process will switch on is a by-chance. We suppose that this 
process is the transition of the neutron star from propeller to 
ejector. Another alternative is the transition closed-open cavern 
around the ejecting neutron star (Zamanov 1995b, Harrison et 
al. 2000). In this case the expected magneto-dipole luminosity 
of the neutron star is about 6 x 10^^ erg s~^, i.e. 10 times more 
than the expected in the ejector-propeller model. 

It can not be excluded that a magnetized black hole is acting 
in LS I+61°303 (Punsly, 1999). If such a compact object does 
exist in LS I+61°303, the microquasar scenario proposed by 
Massi et al. (2001) should be also considered seriously. Further 
observations are required in order to finally solve the true na- 
ture of this X-ray binary. 

7. Conclusions 

We tested the ejector-propeller model for the radio-emitting 
Be/X-ray binary LS 1+61° 303. The calculations show that the 
X-ray variability and radio outbursts can be the result of the 
transition from the propeller to the ejector accretion regimes 
every orbital period. The main results are: 

- The parameters expected from the ejector-propeller model 
are appropriate for the radio plasmon and the calculated ra- 
dio light curves are similar to the observed radio outbursts in 
LS I+61°303. 

- In terms of Ejector-Propeller model the X-ray maximum 
is a result of the propeller action during the periastron passage 
of the neutron star. 



- The periastron is expected to correspond with radio phase 
0.5 using the late radio ephemeris - i.e. the time of the X-ray 

maximum. 

- The observed phase shift between the radio and the X- 
ray maxima is in agreement with the ejector-propeller model of 
LS 1+61 °303. We expect that the X-ray and the radio maxima 
will always peak at different orbital phase. 

- The total X-ray luminosity expected in terms of Ejector 
-Propeller model 1 — 200 x 10"^"^ erg s^^ is in good agreement 
with the observed by ASCA (Paredes et al. 1997) variabihty 
1 - 6 X 10^^ erg s'^ in the 2-10 keV band. 

There are a few systems where it is expected that the neu- 
tron star survives transitions from propeller to accretor and 
backward - A 0538-66, X 0331+53 (Raguzova & Lipunov, 
1998). In our opinion, LS 1+61°303 is probably the first exam- 
ple where we observe the transition from ejector to propeller 
and backward to occur every orbital period. 

During the lives of neutron stars in binary systems, it is be- 
heved that they evolve through different stages - ejector, pro- 
peller, accretor (Lipunov et al. 1994). Studies over systems 
such as LS I+61°303, A 0538-66 and X 0331+53, which prob- 
ably exhibit transitions between these different regimes, will 
give us the unique possibility to investigate these transitions 
over a time scale of months. This is an interval of time ex- 
tremely short considering that the transitions, which are result 
of the neutron star evolution, occur on much longer time scales 
(^10 Myr). New Ught over the evolution of the neutron stars 
in binary systems is expected from this kind of work. 
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(PB97-0903) and by Junta de Andaluda (Spain). 
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